Witnessing Continuous Variable Entanglement with Discrete Measurements 
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In this Letter, we derive an entropic Einstein-Podolsky-Rosen (EPR) steering inequality for con- 
tinuous variable (CV) systems using only experimentally measured discrete probability distributions 
and details of the measurement apparatus. We use this inequality to witness entanglement between 
the positions and momenta of photon pairs generated in spontaneous parametric downconversion 
(SPDC). We examine the asymmetry between parties in this inequality, and show that this asym- 
metry can be used to reduce the technical requirements of experimental setups intended to witness 
entanglement. Furthermore, we develop a more stringent steering inequality that is symmetric 
between parties, and use it to witness symmetric EPR steering. 



Witnessing entanglement is a task whose difficulty 
grows with the dimensionality of the system. Witness- 
ing continuous variable (CV) entanglement in the labora- 
tory is understandably challenging because CV systems 
are infinite-dimensional. As interest in high dimensional 
entanglement increases, including interest regarding its 
fundamental aspects [H-flU- as well as its use in applica- 
tions such as quantum key distribution Q , quantum tele- 
portation and quantum computation witnesses 
of CV entanglement become increasingly important. 

Current experimental witnesses of CV entanglement 
that don't require a complete reconstruction of the den- 
sity operator or rely on prior knowledge of the system 
are based either on coarse-grained approximations to the 
density operator or conditional variance products Q, 
both of which encounter difficulties at larger dimensions. 
Since there are as yet no Bell inequalities [13, [3] for 
CV systems, the strongest level of CV entanglement that 
can currently be witnessed is EPR steering [15| with the 
violation of an EPR steering inequality. 

EPR steering, first formulated by Wiseman et al. 
requires a level of entanglement weaker than Bell nonlo- 
cality [13], but stronger than mere nonseparability [l6| . 
Such steerable states have correlations strong enough to 
rule out any model of local hidden states (LHS) where 
conditional distributions of measurement outcomes can 
be modeled as those arising from single quantum states 
17| . All current steering inequalities arise from demon- 
strations of the EPR paradox [ISJ , when the conditional 
probability distribution of measurement outcomes of a 
given particle B, conditioned on measurements of another 
particle A, is localized to a precision better than the un- 
certainty principle on particle B allows. Since violating 
the single-particle uncertainty principle is forbidden for 
single quantum states, demonstrating the EPR paradox 
necessarily witnesses EPR steering, which is useful since 
it allows any formulation of the uncertainty principle to 
be adapted into an EPR steering inequality 17 [. 

Reid [3] was the first to create an EPR steering inequal- 
ity using the Heisenberg uncertainty principle. Later, 



Walborn et al. [^ developed a steering inequality using 
an entropic [l9| formulation of the uncertainty princi- 
ple for CV position and momentum [i^. For all states 
whose spatial correlations are insufficient to demonstrate 
the EPR paradox, the inequality. 



h{xB\xA) + h{kB\kA) > nlog(7re). 



(1) 



must be satisfied where: n is the number of spatial di- 
mensions; h{xB\xA) is the continuous Shannon entropy 



21| of the distribution of measurement outcomes of xb 
conditioned on outcomes of a;^; h{kB\kA) is similarly de- 
fined for measurements of the wavenumber or momentum 
in natural units; and the base of the logarithm is deter- 
mined by the units in which we choose to measure the 
entropy. Walborn et al. proved this inequality for one 
dimension, but it is easily generalized to n dimensions 
assuming that different spatial degrees of freedom are 
statistically independent. Since the entropic uncertainty 
principle is more inclusive than Heisenberg's uncertainty 
principle the EPR steering inequality derived by 

Walborn et al. is more inclusive than Reid's inequality, 
as Walborn et al. demonstrate with Hermite-Gauss bipar- 
tite wavefunctions Insightful as the entropic steering 
inequality ([T]) is, it cannot be used in the laboratory be- 
cause continuous probability densities cannot be deter- 
mined with a finite number of measurements. 

In this Letter, we derive an exact entropic EPR steer- 
ing inequality suitable for experimental investigations of 
CV position-momentum entanglement with discrete mea- 
surements: 



HIXbIXa) 



H{Kb\Ka)>Y.^oi 



AxBiAkBi 



(2) 



Here we have that H{Xb\Xa) is the discrete Shannon 
entropy [2l| of measurements of the position of particle B 
conditioned on measurements of the position of particle 
A, where both position domains have been discretized 
into equally spaced windows reflecting the precision of 
the experimental setup. 



2 



Our inequality ^ provides a simpler method to suc- 
cessfully witness EPR steering. We do not need to 
reconstruct CV probability density functions, and we 
only require experimentally resolvable discrete probabil- 
ities along with those details of the experimental setup 
needed to determine the measurement resolutions Axsi 
and AkBi- Though Reid's inequality is also experimen- 
tally useful for variances conditioned on particular mea- 
surements, averaging these as a witness for EPR steering 
may fail to violate their inequality unless all possible con- 
ditional variances violate the uncertainty bound. Since 
our inequality examines the joint probability distribution 
of measurement outcomes, we cover all possible condi- 
tional measurements in a comprehensive way. 

To derive our inequality, we will use a fundamental 
connection between continuous and discrete entropies ([8]) 
studied in Ref. 22], to show that any two continuous 
random variables x and y that can be discretized into 
equally spaced windows of size Aa; and Ay satisfy the 
following inequality; 



h{y\x)<H{Y\X)+\og{Ay). 



(3) 



Consider an experiment to measure random variable x 
which can take the value of any real number with prob- 
ability density p{x). The experiment is only capable of 
measuring x to discrete windows Xi of size Aa;. The 
probability of measuring x to be in window Xi is 



P{X,) = / dx p{x) 



(4) 



where the region of integration Aa;^ is the range of values 
of X between Xi — '^Ax and xg+^Ax, and xg is the value of 
X at the center of the window Xg. The Shannon entropy 
of this discrete probability distribution is given by 



HiX) = -J2PiXi)log{P{Xg)), 



(5) 



and the Shannon entropy [21] of the continuous proba- 
bility density function p(x) is expressed as 



h{x) 



dx p{x)log{p{x)). 



(6) 



We now define the distribution pe{x) as the probability 
distribution of x conditioned on being measured within 
window Xg. The continuous entropy hg{x) is defined as 
the entropy of pg{x) where 



pg(x) 



(7) 



for all values of x in the window Xg, and is zero other- 
wise. By breaking up the continuous entropy h{x) into 
a sum over all windows, and expressing h{x) in terms of 



hg{x) and P{Xg), we obtain the fundamental connection 
between discrete and continuous entropies; 



h{x) =J2Pi^Mx) + H{X). 



(8) 



This connection exists for joint entropies as well as 
marginal entropies, where we now define hgm{x, y) as the 
entropy of the joint distribution pgm{x, y) conditioned on 
X being measured within window Xg and y being mea- 
sured within window Ym. 

The conditional entropies h{y\x) and H{Y\X) are de- 
fined as differences between joint and marginal entropies 

in, 



(9a) 
(9b) 



h{y\x) = h{x,y) - h{x), 
H{Y\X) = H{X,Y) - H{X). 



By using ^ for both single and joint entropies and ([9aj) 
and the knowledge that conditioning on additional events 
reduces the average entropy, it can be shown that 



h{y\x) < P{X^, Y^)hg^{y\x) + H{Y\X) 



(10) 



The uniform distribution maximizes the entropy |21| . so 
that when all windows Ay^n are of equal size, we have 
hgm{y\x) < \og{Ay), which completes our proof of ([3]). 
Where XAi and XBi are another ordinary pair of ran- 
dom variables, we can substitute the expression into 
the steering inequality created by Walborn et al. ([T]) to 
derive our entropic EPR steering inequality suitable for 
experimental investigations of CV entanglement. For a 
particular spatial degree of freedom i € {1, (i.e. a 

particular dimension in space) we've shown that 



H{Xb^\Xm) + H{KB^\KA^) > log 



AxBiAksi 



(11) 



When different spatial degrees of freedom are statistically 
independent of one another, the entropies add, giving us 
the n-dimensional discrete steering inequality we sought 
to prove ([2]). 

When applying our steering inequality ([2]), there are 
a number of critical details to consider. First, our dis- 
crete steering inequalities for each degree of freedom pT|) 
have a cutoff of experimental resolution below which our 
ability to probe quantum phenomena dependent on the 
uncertainty principle ceases to exist. For any AxBiAkBi 
larger than 7re, we have insufficient resolution to witness 
entanglement in the i^^ degree of freedom with our steer- 
ing inequality (jlip . In this regime, the steering inequality 
for the i^^ degree of freedom pT|) becomes impossible to 
violate since the bound on the right hand side becomes 
negative, while the sum of discrete entropies on the left 
hand side is always nonnegative. However, as we increase 
the resolution, decreasing AxBiAkBi, we are better able 
to violate our EPR steering inequality. 
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Secondly, the inequality ([2]) appears to depend on the 
resolutions of only one detector. Though the conditional 
entropies have an inherent dependence on the resolution 
of both parties, the inequality is symmetric, as we show 
in Fig. [1] with experimental data from Ref. [5| . 



o 

01 

Dl 
o; 
D 



15 



10 



-5 ■ 



-10 





B = 24x24 




B= 12x12 




^ — ^ = 8jL^ 
— - — — ^^J^ 




-® B = 4x4^ 




'0 

^ — ^^iJ^i 




B = 2x2 



10 15 20 

Pixels per side of DIVID A 



25 



FIG. 1. Violation of our EPR steering inequality ((2]) for mea- 
surements of party B conditioned on measurements of party A 
as measured in standard deviations of the difference between 
the bound and the sum of conditional entropies. This plot 
was formed by independently downsampling the joint distri- 
butions for a 24 X 24 data set. A positive value represents 
a violation. We see that below 4x4 resolution for party B 
as suggested by theory, the inequality is not violated for any 
resolution of DMD A, but violation is possible for large reso- 
lution of party B and party A having resolution below 4x4. 



Using Bayes' rule to swap parties [2l|, we see that sym- 
metry between parties exists only when the marginal en- 
tropies for each party are equal to one another. Due to 
the asymmetry between parties in it is not true in 
general that the ability to violate the discrete inequality 
on a given subsystem is simply limited by the party with 
lowest resolution. However, since the double gaussian 
wavefunction used to model SPDC is symmetric between 
parties experimental investigations of EPR steering 
in SPDC do not exhibit a high degree of asymmetry, 
though the asymmetry is still significant enough to re- 
duce the necessary technical specifications of an exper- 
iment to witness EPR steering. Particularly at higher 
resolutions, one detector can have less than half the res- 
olution of the other and still be used to witness steering. 

In addition, we consider the case where the product 
of position and momentum windows is fixed, but each is 
allowed to vary. For large position window sizes: the po- 
sition conditional entropy decreases toward zero, and the 
momentum conditional entropy increases without limit. 
This makes the sum of discrete entropies in ^ grow with- 
out limit while the uncertainty bound remains constant. 



making EPR steering progressively more difficult to ob- 
serve. Downsampling the momentum distribution would 
not make the inequality easier to violate because doing 
so would also decrease the bound on the right hand side 
even though it would result in a smaller discrete entropy 
on the left hand side of the inequality. 

Expressing our steering inequality in terms of condi- 
tional entropies is useful, but there is also a stronger level 
of entanglement (which we call symmetric EPR steering) 
sufficient to allow EPR steering between both parties. 
This level of entanglement can be witnessed by express- 
ing our inequality in terms of the mutual information. 
The mutual information ^27] is defined as 

I{Xa,Xb) ^ H{Xa) + H{Xb) - H{Xa,Xb). (12) 

Our steering inequality ^ expressed in terms of mutual 
information becomes 



I{Xa. Xb) + I{Ka. Kb) < J2 l^S 

1=1 

+{H{Xb) + H{Kb))- 



AxBiAkBi 



(13) 



Since the entropies under discussion are discrete, they 
are bounded above by the logarithm of the number of 
windows in the viewing area which can be expressed in 
terms of ratios ^i^Lki — where L^i and Lt,- are defined 

Ax Bi^KBi ' 

as the total extent of the viewing area in the i*'* direction 
for position and momentum measurements of party B, 
respectively. Canceling out like terms and simplifying, 
we arrive at a more restrictive steering inequality which 
is symmetric between parties; 



IiXA,XB) + IiKA,KB) <l0g 



(7re)" 



(14) 



Since equation ((Tl)) can also be derived by exchanging 
parties A and B, L^i and L^i can refer to the viewing 
areas of either party. Violation of this symmetric steer- 
ing inequality simultaneously witnesses EPR steering for 
both parties; no conditional distribution of measurement 
outcomes can be ascribed to a single quantum state. 

To use these inequalities in practice, we used data from 
the experiment in Ref. |5| where we performed measure- 
ments of the near field and far field probability distri- 
butions (positions and momenta) of pairs of entangled 
photons generated in SPDC by assembling histograms of 
coincidence counts. In Ref. Q, we were able to use this 
joint probability distribution to calculate the conditional 
entropy and the mutual information, but not to deter- 
mine whether the system had entanglement. With our 
new EPR steering inequalities, we can verify that the sys- 
tem in Ref. [1] at certain resolutions does exhibit both 
EPR and symmetric EPR steering, and therefore is not 
only entangled, but sufficiently entangled to demonstrate 
the EPR paradox. 
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In order to better explain what was measured, we 
provide a brief description of the experimental setup in 
Ref. [5] which generated the joint probability distribution 
that we now analyze with our discrete steering inequal- 
ities (P)) p^ . We separated the downconverted photons 
with a 50:50 beamsplitter into signal and idler arms, and 
measured coincident detections using time correlated sin- 
gle photon counting. With these coincidence counts, we 
measured the joint transverse spatial probability distri- 
butions of the photon pairs by imaging the face of the 
nonlinear crystal (and then its Fourier transform) onto 
DMD (Digital Micromirror Device) arrays in each arm 
which allowed us to look at both spatial degrees of free- 
dom in the transverse plane. Uncertainties in the prob- 
ability distribution were estimated by assuming Poisso- 
nian statistics. We measured the joint probability distri- 
butions at a variety of different resolutions between 8x8 
pixels and 24 x 24 pixels with the same total viewing area. 
The time it took to gather a full set of coincidence counts 
was the primary factor limiting our ability to measure at 
higher resolutions. At larger resolutions, the number of 
joint pixel configurations needed to cover the joint prob- 
ability distribution increases as well as the time it takes 
to measure a single one of these pixel pairs. Primary 
sources of error were due to the temperature instability 
of the nonlinear crystal over the time scales needed to 
take data, the imperfect alignment of the DMD arrays, 
and the imperfect in-coupling of light from the DMD ar- 
rays into our photodetectors. 

With the joint probability distributions obtained from 
measurements for both the positions and momenta of the 
photon pairs, we calculated the conditional entropies that 
go into our steering inequality, and used the details of 
the experimental setup to determine all the window sizes 
Axsi and Afc^i. Given the details of the experimental 
setup, and that the two transverse degrees of freedom 
are considered independent of one another, the discrete 
steering inequality ([2]) takes the form 



H{Xb\Xa)+H{Kb\Ka) > 2 log 



/ 7re 



(15) 



where Aks is the resolution of the detector of party B 
to distinguish differences in the horizontal component of 
the momentum. 

The experimental data from Ref. Q relevant to this in- 
equality is displayed in Fig. [5] In our setup, we were able 
to violate this inequality by between 3.6 and 16.4 stan- 
dard deviations for 8 x 8 resolution to 24 x 24 resolution 
respectively. The values of and Lk used to determine 
Aks and Axb were 1.04 x lO^^m and 1.00 x lO^m"-^, 
respectively. The symmetric steering inequality (|14l) was 
not violated for 8 x 8 resolution, but was violated by 
between 3.4 and 7.0 standard deviations for 16 x 16 res- 
olution and by between 6.6 and 10.7 standard deviations 
for 24 X 24 resolution. 
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FIG. 2. Experimentally determined sums of conditional en- 
tropies in our EPR steering inequality (|15|) as a function of the 
product of resolutions and . The different symbols des- 
ignate different resolutions of experimental data which were 
then downsampled. For example, the triangles represent data 
from an experiment recording at 24 x 24 resolution. Since 
the product of the window sizes is dimensionless, this can be 
expressed in any units we wish. The region below the black 
curve is the region where our EPR steering inequality (|15p is 
violated; the sum is less than the bound. 



In this Letter, we created two new EPR steering in- 
equalities (([U and (IT4)) ) especially well-suited for exper- 
imental investigations of CV position-momentum entan- 
glement. Our first inequality is more inclusive and with 
fewer complications than Reid's inequality Q at suffi- 
ciently high resolution, while both are easier to imple- 
ment in the lab than the inequality derived by Walborn 
et al. H)). We have successfully witnessed both EPR 
steering and symmetric EPR steering with experimental 
data from Ref. 5], and showed that there is a demon- 
strable asymmetry between parties which allows steer- 
ing to be witnessed even when one detector has compar- 
atively low resolution provided the other is sufficiently 
high. These inequalities are powerful and effective char- 
acterization tools that we expect to be widely used in ap- 
plications ranging from future quantum communication 
networks to fundamental physical experiments involving 
high-dimensional quantum states. 
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